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Abstract 

In present study I examine the capability of diffusive shock acceleration mechanism to explain 
existing data on radio emission from evolved large diameter shell-type adiabatic supernova remnants 
(SNRs). Time-dependent "onion-shell" model for the radio emission of SNRs is developed, which 
is based on the assumptions: a) acceleration takes place from thermal energies and test-particle 
approximation is valid; b) the problem of injection is avoided by introducing, like Bell (1978), two 
injection parameters; c) to take into consideration very late stages of SNR evolution the analytic 
approximation of Cox and Andersen (1982) for the shell structure is used; c) no radiative cooling. 
Constructed Surface Brightness - Diameter (E — D) tracks are compared with the empirical £ — D 
diagram. The main conclusion of the study is that the DSA mechanism is capable of explaining all 
the statistics of radio SNRs including very large diameter remnants and giant galactic loops 

1 Introduction: 

Diffusive shock acceleration (DSA) mechanism are believed to be the source of radio emitting 
relativistic electrons in young and middle age SNRs of diameter D < 20 pc. But for the radio 
emission of evolved large diameter SNRs the mechanism of van der Laan (1962) or some modification 
of this mechanism (e.g., Blandford and Cowie 1982) are widely believed to be responsible. According 
to these models preexisting in the ambient ISM electrons compressed to high level at the radiative 
shock are responsible for the radio emission of the remnant. Due to the unstability of radiative 
shock waves this mechanism may however encounter a number of difficulties in reproducing the radio 
emission from very large diameter SNRs evolving in warm phase of the ISM. Moreover, if the density 
of the ambient medium is very low, the SNR will finish its life by merging with the ISM before cooling 
becomes important. For such remnants DSA becomes a main candidate for generation of relativistic 
electrons radio emitting in the magnetic fields of 10 -5 -f- 10~ 6 G. But electron acceleration in SNRs 
is more difficult to estimate quantitatively since the injection efficiency and even the very process of 
acceleration for electrons are still unclear. 

The goal of present study is to apply DSA mechanism, under very simple and common assumptions 
about the injection, to follow the evolution of shell-type SNRs up to the very large radii and, at the 
same time, to obtain some new constraints on the theory of acceleration mechanism if we will be able 
to achieve agreement between the model predictions and observations. 

2 The Model: 

We used the onion-shell model of Moraal and Axford (1983) as it was done in Asvarov (1992, 
1994) where the Sedov solution for the remnant structure had been adopted. In present study SNR 
is modeled by using an analytical approximation of Cox and Anderson (1982) which follows the 
development of an adiabatic spherical blast wave in homogeneous ambient medium of finite pressure. 
At early times this approximation resembles the zero pressure Sedov similarity solution but extends 
the range of investigation well into the regime in which the external pressure is significant. 

To avoid the problem of injection we introduced, like Bell (1978), two injection parameters: we 
inject electrons as "test particles" at momentum pi n j = -0Pth ; where p t h = (2m e T s ) 1//2 is the down- 



stream electron thermal momentum, the concentration 01 which assumed to be proportional to the 
density of the ambient thermal electrons: iVi n j = ipn oc . At the shock front equipartition between 
electron and proton temperatures is assumed which means that p- m j is proportional to the velocity of 
shock wave. 

For the accelerated electrons from the energy loss processes only adiabatic cooling was taken into 
consideration. 

Assuming the magnetic field to be frozen in to the plasma, we model the density dependence of 
H as H = Hq(p/ '/9o) k > where Hq is the ambient value of the magnetic field strength. 

Although in our analysis we consider only extended SNRs we calculate several models in high 
density environments for which the radiative phase begins relatively soon after the explosion of SN. 
In this case total flux is obtained by integration over the layers of the shell where radiative cooling 
does not occur. This implies that we completely ignore the action of DSA at radiative shocks. 

3 Results and Discussion: 

Empirical T, — D relations are very useful tools for testing the theoretical models of SNR evolution. 
Before comparing our model with the observations, we formulate the common properties of the model 
predictions. Here we concentrate on two main observable radio characteristics of SNRs: the spectral 
index and surface brightness. DSA at strong shock waves in test particle approximation predicts for 
spectral index the value of 0.5, which will increase as the shock intensity (Mach number) decreases. 
Calculations show that at Mach numbers M < 4 the value of the mean radio spectral index gets 
greater then 0.6 but it remains bounded by maximal value of 0.75 during the following evolution. 
This is the result of distribution of magnetic field strength at the shock but the real average spectrum 
of electrons inside the remnant will be somewhat softer than 2a + 1 = 2.5. 

What concerns another important radio characteristic of the SNR, the surface brightness, the rem- 
nant evolves at nearly constant radio surface brightness followed by very steep drop. It is important 
to note that the dependence of these and other radio characteristics of the remnant on shock Mach 
number has almost universal nature, very weakly depending on the input parameters. 

To compare the model predictions with the observations we have collected a set of shell-type and 
a few composite SNRs in our Galaxy with known distances and remnants in LMC. Corresponding 
£ — D diagram is shown in Fig.l. Error bars are due to uncertainties in the distances for some SNRs; 
for several objects only lower and for one SNR upper limits are known. 

As a standard set of input parameters we used: the energy of SN explosion Est? which is varied 
in the range (1 -j- 5) 10 51 erg; the ambient thermal electrons density nQ e in the range (5 -j- 5 10~ 3 ) 
cm~ 3 ; the strength of the ambient magnetic field Hq in the range(3 -j- 10) x 10 _6 G. In all models the 
mass of SN ejecta is taken to be one solar mass. In Fig.l several Eighz(-D) - tracks are shown for 
different values of n oe , -E?sn and Hq. In all calculations we used tp = 3 and ip = 4 x 10~ 4 justifying 
the use of test particle approximation. It is important to note that the shapes of the evolutionary 
tracks depend very weakly on these parameters although the magnitude of £ depends on (p linearly. 

As can be seen in Fig.l varying mainly the values of n oc and E^n the model is able to cover all the 
remnants in the £ — D diagram including two large diameter relatively bright SNRs (HB 9, OA 184) 
and giant radio loops. As the latter is attained with the price of somewhat large value for i?sN of 
5 x 10 51 erg we calculated several models in which the contribution of the ambient relativistic electrons 
was included via the injection of them into the DSA. As an example we have taken the spectrum 
j = 1.4 x W 2 E~ 2 - 2 electrons m" 2 s" 1 sr" 1 GeV" 1 (E in GeV) from Fichtel et al. (1991). In Fig.l 
two tracks are drawn by dashed lines. 

It is interesting to not that the shapes of predicted by our model evolutionary tracks are in 
excellent accordance with the prediction made by Berkhuijsen (1986) that "radio remnants may 
evolve adiabatically at nearly constant Sr, followed by a steep decrease" 




Figure 1: Surface brightness - diameter diagram at 1 GHz for SNRs in Galaxy (circles) and LMC 
(squares) and giant galactic loops (crosses). All the data for galactic SNRs are taken from the 
catalogue of Green (1998), the data for Loops are taken from Berkhuijsen (1986), and the data for 
SNRs in LMC are from (Mills et al. 1984). Galactic SNRs with a < 0.4 or/and composite SNRs 
indicated by open circles. Curves are the modeled evolutionary tracks: 1 - n oc =5, Hq=1, Esn=2; 2 
- the same as 1 but E$n = 1; 3 - n oc = 1, Hq = 1, £^sn = 1; 4 - n oc = 0.05, Hq = 0.5, £"sn = 2; 5 - 
n oe = 0.01, H = 0.3, £ S N = 5; 6 - n oc = 5 10" 3 , Hq = 0.3, E SN = 2; 7 - n oc = 5 10~ 3 , Hq = 0.3, 
Ssn = 1- Here -Esn is in 10 51 erg, Hq is in 10~ 5 G, n oc is in cm -3 . Dashed curves correspond to 
models in which the contribution of background cosmic ray electrons was considered. 

As can be seen in Fig.l radio evolution of SNRs depends on two parameters, £sn and no e , equally. 

In the framework of present model a number of features of the empirical S — D relation obtains 
a simple explanation. For instance, the small number of remnants with small diameters and low 
£ (lower left corner in the diagram) is the result of very fast evolution of SN blast wave in the 
low density ISM where SNRs have low S. According to our model it is easy to account for high 
concentration of SNRs at diameters 30 — 50 pc in the £ — D diagram: different kind of evolutionary 
tracks intersect at these diameters and the sample of remnants here consists of objects evolving at 
different initial conditions. Of course, in the origin of the empirical £ — D relations we can not exclude 
at all the contribution of various selection effects. Moreover, not all the remnants can be described 
by our model. Indeed, in the £ — D diagram (Fig.l) the composite SNRs and SNRs with a < 0.4 
(indicated by open circles) have systematically large values of £ , which can be understood that in 
these remnants an additional more effective mechanism acts. 

It is important to note that adopted values for Bells' parameters do not contradict the observations 
at standard values for the input parameters, characterizing the ISM and the SNR itself. This fact 
implies that the test particle approximation has factual realization in evolved shell- type SNRs. One 
more argument which favors the DSA mechanism is the statistics of spectral indices. The catalogue 
of SNRs of Green (1998) contains 80 SNRs with well determined values of the spectral indices,a, 
from which 57 remnants (71%) have a > 0.45 and only two SNRs (one of them is young peculiar 
SNRs Cas A) a > 0.75. Practically there are no objects contradicting the prediction of our model 
that a max < 0.75. It is well known that young SNRs have systematically large values of a, which can 
be explained by the back reaction effects or by the action of other then DSA mechanisms. According 



to the model blNKs with Mach numbers M < 4 has the mean value ot a > 0.6. Assuming tor 
simplicity that SNR evolves according to the Sedov law, M oc i _3//5 , for the number of SNRs with 
Mach numbers greater then M we have N(> M) oc M -5 / 3 from which it follows that the number of 
SNRs with a > 0.6, N(a > 0.6) = N(M < 4), must be about 55% of total number of SNRs. Here 
we have adopted for the final Mach number Mf = 2.5 at which £ drops more then two order of its 
initial value and the SNR becomes invisible. In the catalogue Green (1998) 22 out of 80 SNRs have 
0.60 < a < 0.75 which makes 27.5%. This discrepancy easily can be explained by the selection effects 
that SNRs with large diameters and small Mach numbers have low £, consequently, difficult to be 
detected, though their number is more than the number of bright remnants by a factor of 8 — 10. 

The mean size and age of SNR with E$n = 10 51 ergs, evolving in the ISM with nQ e = 0.01 cm~ 3 , 
when S drops to ~ 5 x 10" 22 (Wm- 2 sr- 1 Hz- 1 ) are 150 pc and 2 x 10 years, respectively. If 
SN occur with a rate of 30 year -1 then such SNRs occupy 3.5 x 10 65 cm 3 of the galactic volume 
of 7T x (25kpc) 2 x (lkpc)/6 = 9.65 x 10 66 cm 3 , or 1 part in 28(3.6%). This estimate depends on 
the value of £sn as E^. The probability that random line of sight will hit such SNR is 0.21 or 
10 in 47 and the dependence on E$n is the same. The last estimate shows that SNRs can play 
important role in the origin of background radio and gamma emissions of our Galaxy. Predicted by 
our model spectral indices are in accordance with the radio background observations. What concerns 
the gamma-background our model in accordance with the model of "bubbling swiss cheese" of Pohl 
& Esposito (1998) thought the value of electron spectral index of 2.0 demanded in their model is in 
contradiction with our predictions. 

4 Conclusions: 

The model based on the assumption that the radio emitting electrons are accelerated by DSAmech- 
anism very well explains the statistics of shell-type radio SNRs. From this we can conclude that the 
test particle approximation and the supposition that acceleration of electrons takes place from the 
thermal energies has realization in evolved SNRs in spite of all theoreticaldifficulties concerning the 
physics of collisionless shock waves. 

We obtain that ip = Pmj/pth — 3. The idea that there is no acceleration at the radiative stage of 
SNR evolution does not contradict observation. 

Presented model also in accordance with radio and gamma background observations. 
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